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ABSTRACT 
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A cylindrical solar water heater has been fabricated, 
theoretically analysed and experimented upon in the present 
study. The heater consists of two glass tubes of 15.0 and 
17.8 cei diameter enclosing a G-.I. pipe of 11 era'*, diameter# 

Water is made to flow in the G.I. tube whose outer surface is 
painted black, which acts as a solar energy absorber. 

Experiments have been carried cut to test the perfor- 
mance of the water heater under two different modes of operation, 
namely 

a) natural oonvection 

b) forced flow condition for mass flow rates of 
30, 90, 150 and 210 kg/h. 

The day-long collection efficiency under the first mode 
has been found to he 31.5$ while in the second mode it lies bet- 
ween 62 - 77$. The maximum efficiency of 76.4$ has been attained 
for mass flow rate of 150 kg/h. 

An approximate theoretical analysis has also been carried 
out on the heater for forced flow conditions. The experimental 
results are compared with -those obtained "theoretically, and the 
comparison shows a reasonably good degree of agreement. 



CHAPTER I 


INTROHJCTION 

1.1 HEED KK SOLAR ENERGY 
UTILIZATION 

With, the increase in industrialization and advancement 
of technology there has been an evergrowing demand for energy 
both in the developed and the underdeveloped countries. To get 
some idea of the rate of consumption of energy we shall quantify 
power from 33,000 million tons of soft coal by 1 Q * [_lT . Until 
1850, man had used 9Q. In the next 50 years man had used another 
4Q. The present rate of consumption averages to about 10Q per 
century. By 1'980 we shall be using IQ every five years. This 
means roughly speaking, that half of all the energy consumed by 
man in the past 2000 years has been consumed in the last me 
hundred. 

Together with this catastrophe of dwindling fossil 
fuels, is the tempering of ecological balance of nature because 
of pollution which comes in the wake of consumption of fossil 
fuels. Thus there is need to look fbr alternative sources of 
inexhaustible and non polluting energy. 

Among the alternative energy resources that are being 
seriously,; considered at present are j i) Energy from fast 
breeder reactors, ii) Geothermal energy and iii) Solar Energy. 

It is believed £2j that if the fast breeder reactor technology 
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comes through, then the existing -uranium resources of the world 
may last another 1,000 years from now. But this has the serious 
disadvantage of ecological pollution and scientists have still 
not been able to dispose off radioactive waste safely. &eo- 
thermal Energy also offers an alternative to the Ehergy crisis 
but at best the stored thermal energy in the world’s major 
geothermal areas would, if completely exploited, provide energy 
for about 50 years. Even in that time the contribution to the 
world’ s annual consumpt ion would be le ss than that from tidal 
power which if fully exploited is about 1 $ of the world's 
Energy consumption of 1970. 

Thus we come to the 3rd alternative i.e. the Solar 
Energy. It has its own drawbacks like capital cost of equip- 
ment etc. but from availability point of view it is inexhaustible 
and what is more - non polluting. Thus it becomes quite obvious 
that this Energy be tapped and used for future power needs. 

1.2 AVAILABILITY OF SO LAB ENERGY 

The amount of Solar Energy incident outside the atmos- 

1A ■ r ■ ■ 

phere amounts to 1.5 x 10 kwh/year i 3! . But very little 

V ■' 

of this Efaepgy is available for use by man. Nearly a third of 
it is lost by reflection to space. Another large fraction is 
intercepted by the meteorological thermodynamic engine and its 
associated cloud belts aid cloudy regions. Part of it is 
incident on 70$ of earth’s surface covered by water. Not all 
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of what is left of the solar energy is available for human use. 

IS 

Still tile amount of energy falling is about 0.2 x 10 kwh/year 
which is almost 16,700 times the amount of energy than we now 
consume f4 7 . 

^ - 9 * 

There are two obvious problems in efficiently harness- 
ing solar energy. First, it is not available continuously and 
thus it necessitates the use of some storage device. Secondly, 
the energy is diffused in form. Vftiile the amount of energy 
available is enormous, that available at any point is not large 
enough (solar constant being 2 Cal/cm -min.) to be useful. Thus 
in order that this energy be of use to man it must be efficiently 
and economically harnessed. 

1.3 HAHHESSBTG OF SOLAS MERGY 

Primarily Miere are three ways of tapping solar energy, 
depending upon the form in which we want the output. They are t 
i) Direct Energy Conversion; ii) Concentration of Solar Energy and 
iii) Collection of Solar Energy. The first one is used for 
directly getting electric power. Its primary use has been in 
space stations and satellites. The major drawback in such devices 
has been their low efficiencies like 10 - 15$ and very high cost 
(on the order of Rs. 4 lakhs/kw) . The second method, namely the 
concentration of solar energy is used for getting high temperar- 
tures e.g. for generating steam to run pumps or generate power 
or to use it for studying high temperature technology. Fresenel 
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lenses and parabolic concentrators are normally used for this 
purpose. The major drawback in these types of concentrators is 
the difficulty in making them, their high cost and the rather 
complicated system needed to make them track -the sun. The third 
method, viz. the collection of solar energy is the cheapest aid 
the easiest way of tapping this energy, particularly where the 
requirement of temperatures is of the order of 100 °C. Apart 
from their usage in space heating or cooling, water distillation 
and fruit drying, their maximum use has been for solar water 
heating. This can be used in western countries and in some 
northern regions of India in winter. Apart from this usage, the 
development of these collectors is useful for a country like 
India from the point of view of making an efficient device for 
generating steam or evaporating a liquid of low boiling point 
which in turn will drive a pump. Recently Mienel and Mien el f_5J 
have found that it is feasible as well as 1 economical’ to use 
flat plate selectively absorbing collectors to generate power on 

a large scale. They envisage the total collecting area of 118&.; 

2 6 

118 km which will generate 10 megawatts. Secondly these sort 
of collectors can be used as generators for absorption type of 
refrigeration cycle. Thus to start with there is a need to deve- 
lop efficient solar water heater. 

1.4 REVIEW 03? PREVIOUS WORK 

I Since maximum work bias been done in solar water heat- 
ing, the literature abounds in differ oat types of heatersy/ These 
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heaters differ from one another on minor points but they can 
broadly be classified into three categories. These are : 

i) Tube type h eater : 

This is the most common type of heater. I tb cross 
section is shown in Fig. 1.1 



Fig. 1,1 Fig. 1.2 

Yellot and Sobotka f 6") have investigated -the perfor- 


mance of this type of solar water heater. The steel pipes are 

simply pressed into thermal contact with the absorbing plate. 

The absorbing plate was coated with a selectively absorbing 

coating as developed by Tabor [ 7 ] . There were 2 glass covers. 

2 

The collector area was 2.8 m and steel pipes were of 12 m.m. 
diameter. They have reported an efficiency of 38.5 $ for the 
whole day (natural convection). For flow rates of 44.5 kg/h and 
84 kg/li the average daily collection efficiency was 44.3$ a ^d 
45*5$ respectively. 
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ii) Corrugated sheet type heater 

The heater, Pig. 1.2, consists of a corrugated sheet, 

coated with black paint, attached to a flat sheet by means of 

rivets,. Water flows in the 'pipes' formed by the daeets. Rao 

and Suri £ Sj have experimented with this type of heater. Por 

2 

a collector surface area of 2.32 m and the tank capacity of 136 
litres the rise in temperature of water in winter at Roorkee, in 
natural circulation, v©s from 6 °C to 50 °C (max.). The efficiency 
was found to be 63$. The major problem, in such type of heaters 
is the leakage from the riveted joints. 

iii) Heater cum storage type 

These types of heaters do away with the storage tanks. 

The heater acts as both. Thus the cost is considerably reduced 

in this case. Ghinappa and Gnanalingam ( 9j have made arid tested 

one such heater.; The heater consists of a square coil of 7.55 cm V 

diameter and 13.5 m. long in a wooden box with insulation at J 

2 

bottom and 2 glass covers. Glass surface area is 1.86 m . By 
drawing off water viiai its temperature reaches 49 °C, it is possi- 
ble to obtain 136 - 230 litres of water/day. The collector 
efficiency has been i>und to be 46$. 

Another type of built-in storage water heater has been 
built and tested at Central Arid Zone Research Institute, Jodhpur 
by Garg and Krishnen f 10 1 « 03a e heater consists of a galvanised 

iron (20 gauge) rectangular tank of dimaasions 112 x 80 x 10 can 
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with a capacity of about 90 litres. This tank is placed in a 
rectangular mild sheet tray and insulated from sides and back. 

In front is ordinary glass. The front face of the tank is 
blackaied by lamp black paint. Water temperatures of 55 °C have 
been reported for winter. The efficiency of this type of heater 
has been reported as 75$. 

Chauhan ill j has also built almost the same type of 
heater as above. The capacity of this heater is 70 litres* The 
maximum temperature attained in the month of April has been 
86 °C at about 3*50 p.m. The collector efficiency of 72$ has been 
reported in the forced flow conditions with the mass flow rate of 
76 kg/hr. In the natural convection mode the collector efficiency 
has been found to be between 50 - 58$. 

1.5 MOTIVATION FOE THE PRESENT TOEK 

There are certain disadvantages with the types of solar 
heaters cited above. One of the important disadvantages is the 
flat surface of the collector itself. Set at a specified angle 
of inclination, the surface receives an effective sunbeam, the 
width of which varies between zero and that of the collector 
during the day, as the sun' s position in the sky changes. In 
addition, there exists many avenues for energy losses due to s 
i) conduction from the back and edges of plate, 
requiring heavy insulation for minimization, 
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ii) Radiation and convection from the flat - plate 
surface and the water tubes* 

iii) Reflection from and absorption by glass covers, 
caused by variation in the angle of solar inci- 
dence, and 

iv) induction -through the insulation surrounding the 
pipes connecting the collector and the water tank. 

As a result of these losses, the average yearly collec- 
tor efficiency is between 30 - 40$. In order to improve the 
collector efficiency, a cylindrical solar water heater, using a 
cylindrical surface has been fabricated and investigated by 
Vincze 12 ) . The collector is a tube painted black and covered 

by 2 glass tubes with 2,5 cm space between them. The capacity 

2 

of heater is 22.2 1 and the projected area is 0.186 m . When 
placed at an angle of 63 0 to the horizontal temperatures as 
high as 58 °C have been reached during the' month of August. 

Vincze has reported efficiencies of the order of 216$. The rea- 
son for such high efficiency is given in Chapter IV. The water 
is in the natural circulation. 

As the efficiency appears to be very high it has beei 
thought appropriate to investigate this type of heater more 
thoroughly, especially since it appears to be a very good answer 
for airconditioning devices driven by solar energy. 
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The most important advantage of this type of collector 
is that the effective beam width, for incident radiation remains 
equal to the tube diameter over all the sunshine hours of the day. 
Moreover losses due to conduction at the edges are negligible, 
since the area for axial energy conduction is small compared with 
the tube surface and may be well insulated. 

Thus in the present work a solar water heater of cylin- 
drical type has been fabricated. It consists of a G.I. tube 
of 10.15 can': bore and 11 cm. outer diameter coated with lamp 
black paint. The tube is surrounded by 2 glass tubes of Pyrex 
glass of diameters 15.0 and 17.8 cms. respectively. The capacity 
of the heater is 22.45 litres. The heated section is 1.2 m 
long. 

The heater has been tested for different mass flow 
rates of water starting from 0.5 kg/mi n to 3.5 kg/mi n. It has 
also been tested for natural circulation. The difference between 
the inlet and outlet temperatures of the water has bear noted for 
the whole day and efficiencies have been calculated. The effi- 
ciencies have been found to lie between 62 - 77$ for different 
mass flow rates in ihe range of experimental data. 

An approximate theoretical analysis on this heater has 
also been carried out and the results are compared with those 
obtained experimentally. 



CHAPTER II 


THEORETICAL ANALYSIS 

The cylindrical solar water heater under analysis, 
is as shown in Fig. 2.1. It consists of a metallic tube, which 
has been assumed to be coated with a selectively absorbing 
coating } 7j , placed coaxially between two cylindrical shells 
of glass, which act as energy shields. Each inner tube diame- 
ter is assumed to be J3/2 of that of the tube immediately 
surrounding it, so that most of the radiation transmitted by 
the outer tube may be Incident on the inner tube. The axis of 
the tube is maintained at an appropriate angle of tilt in the 
North - South direction, to effect noimal solar - radiation 
incidence at noon time, on the design day. Water is forced to 
flow at a constant rate through the metallic tube by some exter- 
nal means. 

It is clear that the total incident intensity of ra- 
diation on the tube surface, q, varies with time of the day. as 
a result, the temperature of the metallic tube at any section 
is expected to vary circumferentially as a function of time. In 
order to reduce the complications that will be intro diced in the 
analysis of the temporal and circumferential variations in solar 
energy input, it has been assumed that the total radiation inci- 
dent an the tube surface remains constant at a value equal to the 
average intensity during the sunshine hours, for the place under 
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consideration. This assumption permits us to treat the system 
as one in steady state. Further, the average intensity is 
assumed to be due to a combination of the effects of direct 
solar radiation, diffuse sky radiation and "the reflected ground 
and other long wave length radiations . The equations used in 
this calculation are exhibited later on. 

Since the transmittivity and absorptivity of glass 
vary with the angle of incidence of the radiation, average values 
of these quantities, assumed to remain constant over the whole 
periphery of the glass shields, have to be used in the analysis. 

With the simplifications listed . above, it is possible 
to write down the energy balance equations for the whole system, 

A few other simplifying assumptions used in the analysis are : 

i) The thermal conductivities of the metallic tube 
and glass surfaces are sufficiently high and the 
tube thicknesses so snail, that the thermal resis- 
tances due to the three walls may be neglected; 

ii) Axial conduction through the fluid and the metallic 
tube is negligible, since the product, Re Pr 
(Re = Reynold's number based on the tube diameter and 
■the mean velocity of the flow through the tube, 

Pr = Prandtl number for fluid) is large, being 1280 
for lowest mass flow rate of 0.5 kg/min. and the wall 
-thickness is small; 
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iii) Conduction losses from the ends of the tubes are 
negligible, since the diameter of the outermost 
glass tube is small compared with the solar water 
heater length and the ends are well insulated, 
and 

iv) The fluid properties are invariant in tempera- 
ture because of the small temperature fange 
involved in the heating process. Moreover, the 
viscous dissipation effects have been neglected 
in the formulation of energy balances because 
of Ihe small velocity gradients involved. 

2. 1 ENERGY BALANCE EUQATTONS 

Considering the outer glass tube, it is seen that the 
tube receives energy due to irradiation from the surroundings 
and due to convective as well as radiative energy exchanges with 
the inner glass tube. It looses energy to the environment 
through convection and radiation. 

An energy balance on the outer glass tube therefore 

yields s 

D 1 U i -2 ^ T 1 ~ + D 1 ^ T 1 ” T 2 ^ + D 2 q 

' h a \ (I 2 + (I 2 4 -T a 4 ) (D 

vfoere and are respectively the diameters of the inner and 
outer glass tubes, ^ I* 16 overall heat transfer coefficient 
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between the outer end inner glass tubes and F* „ is the "view 
factor for radiation between the inner and outer glass shields. 
The temperatures of the inner and outer glass shields, respec- 
tively, are T„ and T~, while h , is the heat transfer coefficient 
between the outer cylinder of glass and the ambient air and T„ 
is the ambient temperature. The term ct“ , is the Stefan - 
Boltzman constant. The quantities oC 2 and t 2 are respec- 
tively the absorptivity and emissivity of the outer glass surface 
for solar radiation and long wave length radiation. 

For inner glass tube, similar considerations of energy 
balance yield : 

D 1 U 1-2 (T 1 “ V + <7 *" ^ 1-2 D 1 (T / " T 2 4) 

* 06 1 " C 2 D i * +33 s U s-1 (T s- T 1> + 

+ rr-B F (T 4 - T. 4 ) (2) 

s s— 1 s 1 ' 

where cCj and are respectively the absorptivity and trans- 
mittivity of the inner and outer glass shields, D_ is the metal 

o 

tube diameter. U . is the overall heat transfer coefficient 

T s— 1 

between the metal tube and the inner glass tube and F g __.j is the 
view factor for radiation between the metal tube surface and 
the inner glass tube. 

In writing Eqs, (l) and (2) inter-reflections between 
the inner and outer glass tubes have been considered as small and 
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neglected. Neglect of these interreflections tends to increase 
the apparent energy losses fron the metal tube surface, and 
hence gives a conser-vative estimate for solar heater efficiency* 

Since the liquid flowing into the tube is one of rela- 
tively high Prandtl number (4 Pr 6), it has been assumed 
that the velocity profile is fully developed at the entry of 
the heated section, permitting use of the thermal entry length 
solution for arbitrary wall energy flux J 13 J ■ . For the tube 

surface temperature T we can therefore write the equation : 

s 

x* 

T s &*)' - T e = [d s / (2 k f ) | jf g <x* - t )• Oq -(J)- dt 

0 ( 3 ) 

where, 

2 2 

g (x*> = 4 + exp ( - V ' X* ) / ( V A . ) ( 4 ) 

nn 

and 

q o {x*> = -Cq - U s _ 1 (T s - 3^) - CTF^ (t/ - 30 (5) 

represents the effective energy flux absorbed by tie metal tube 
surface, assuming that the absorptivity is unity. In Eqs. ( 3 ), 
( 4 ) and ( 5 ), 3? is the fluid temperature at the tube inlet and 
x* is a non dimen sionalized axial coordinate measured from the 
tube inlet, given by the equation 

x* = 2 x / (D Re Pr) (6) 

s 

On differentiating eq. ( 3 ) with respect to x* and 
using the rule for differentiation under the integral sign jj14j , 
we obtain, 



16 


r 


X* 


(2 k f / B 3 ) d I g / d x* =J >- 

o 

+ g (o)- q Q ( # } (7) 

If we now substitute for g (x* ~ snd 4 # )- from 

eqs. ( 4 ) and ( 5 ) into the above and use the integral mean value 
theor era. 5 j , we obtain : 

(2 kf / D s ) d T s / d x* = r Ci |x q _ u s _ 1 (l B ( x *) - !_,( x *)) 
~ cr F s-1 [ T s 4 4x V ••(■ x *>j ]• + C,, ^Cq - ^ 

(Tg 4 2C *> - T 1 4 x *>) - <r-I' s _ 1 £ T s 4 <- X* ) - 


where, 


a»d 


1 " S [e*P ( - V n x» ) - 1 ] / (V A n ) 

n 


( 8 ) 

(9) 

( 10 ) 


0 o = 4 + X 1 / (>* 2 A ) 

2 ' £; — ' v n n' 

n 

— * 

while x , represents a point in the interval 
. — * . * 

OCx C_ x , at which the various temperatures are to be eva- 
luated. Further, in eqs. ( 4 ), (9) and (lO), the quantities 
are eigen values and the quantities are constants, listed in 

L«] • 

The mixed mean temperature, T , of the fluid flowing 
through the tube is given by equation; j 13 ( • 

d T q / d x* - (2 D s / k f ) q Q 4 x* >• (ll) 

being the fluid thermal conductivity. 
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If eqs. (2)5 (3), (8) and (ll) are now solved simul- 
taneously to determine the temperatures, T.., T oJ T and T , as 

I c. s in 

functions of x, we can evaluate the heater efficiency (between 
the inlet and the point under consideration) by writing s 

V = a C p ( T q { x ) - tJ / (d 2 x q) (12) 

* 

where m, is the mass flow rate of fluid flowing through the 
tube, and T g is the entry temperature of the fluid. 

2.2 EVALUATION OE CONSTANTS AND HEAT ELUX q 

Before the solution of the above equations be attempted 
it is necessary to find the constants like U.^ gj T J ^ 
etc and the solar radiation input. We shall calculate each of 
these one by one. 

Solar Energy Radiation Input q 

The total radiation falling on a horizontal plane at 
I.I.T. Kanpur has been recorded for every day by the Civil 
Engineering Department on the pyrailometer supplied by Indian 
Meteorological Department. Graphs have been taken from them both 
for direct and diffuse radiation. 

The total radiation on the collector consists of; j 

i) direct radiation on the top surface, ii) diffuse sky radia- 
tion on the top surface and iii) ground reflected radiation on 
the bottom surface. 

‘ TlT KANPUR 
CENTRAL LIBRAE 

Acc. No. A 3&8-I 
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The heater is inclined, to the horizontal at an angle 

r* i 

calculated from the equations : M6 j 

Sin f> = Cos L Cos b Cos H + Sin 1 Sin £ 

Sin (/ = Cos S Sin H / Cos , and (l3) 

Cos 6 = Cos p Cos jzf Sin ^ + Sin Cos 

where 1 = latitude angle; for Kanpur it is 26°26', £= declina- ' 
tion angle and H = house angle = 54* . The sunshine hours are 

taken to be from 7 AM to 5 PM. It has been assumed that the de- 

clination angle changes very little between 29th May to June 12 
(test dates)} ranging between 21° and 23°. Thus the average angle 
has been taken to be 22°. 

Assuming that the direct solar radiation should be normal 
to Ihe tube axis at 12 noon for maximum utilization on the design 
day (i.e. 9 = 0 or Cos 0 = l), the equations above yield; 

VL * jgj go 

which is the inclination of the tube to the horizontal. 

Thus the direct normal radiation, as seen from the 
graph of insolation, be multiplied by l/CosJp to make it normal 
on the inclined surface. 

The diffuse radiation falling on the inclined surface 
should be multiplied by a factor : | 17] 

R = (l + COS If ) / 2 

s * 

so as to make it normal on the inclined surface. 



19 


The ground reflected radiation, falling only on the 
under side of the heater, is found as follows : 


R, 


'GR 


( R 


IN 


+ V 


R 


D g 


■where R^ = Direct normal radiation, R^ = normal diffuse ra- 
diation, = shape factor between the ground and the heater = 

.81 and r = ground reflectance = .12 for gravel roof 

6 

top. Hs] 

The total radiation, for the day of experiment, falling 
on the heater can thus be calculated as : 

i) Direct radiation • This radiation is assumed to 
fall only on the projected area formed by an arc 
of 120° of the cylinder. The rest gets all re- 
flected from the surface. This area for the heater 
2 

is Ov 186 m . 

ii) Diffuse radiation is assumed to fall only on the 

2 

upper portion and the projected area is 0.214 m 

iii) Ground reflected reflection is assumed to fall only 

on the under side of the collector and the projec- 

2 

ted area is 0.214 m . 

Overall Heat Transfer Coefficients 

The overall heat transfer coefficients, U.. „ and U 

1-*! S— 1 

appearing in eqs. (2), ( 3 ) and (8), have to be obtained from the 
curves on overall film conductances for free convection flow 
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between concentric cylinders provided by Jakob | 19 ] . Under 
the conditions of the present computations where Dg = (2 / «/3 ) 

= (4/3) D j the overall film coefficients are very nearly equal 

O 

to the respective thermal resistances caused by the presence of 
annular air gaps, we can therefore write for U and U the 
approximate equations ; 

B 1 D 1-2 = 2 \ / to (D 2 / I,), M D s n s _ 1 = !i a /b(j,/ 3 S ) 

(15) 

Outside Heat Transfer Coefficient 


Assuming the average wind velocity of 12*15 km/h in the 
East-West direction for Kanpur j 20! the heat transfer coeffi- 
cient h Q , between the outside air and outer glass shield of dia- 
meter, Dg can be calculated by using the equation'. \_2 ll 


h D 
a 2 


- a ( Re y 


where k is the thermal conductivity of the ambient air 
8 . 

and the Reynold's number, Re is based on the outer tube diameter 


and the design wind velocity. The quantities, a and b, are cons- 


tants which depend upon Reynold's number. The value of h & , cal- 

2 

culated in this manner comes out to be 15 kcal/hr - m - °C. 


Shape Pactors 

The calculations involve shape factors for a system of 
two concentric and coaxial cylinders, where the inner surface of 
the outer cylinder radiates to itself. Shape factors for this 
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case has been presented in the form of graphs by Hottel and 
Sarofim j_22] and have been used in the present calculations. 

A knowledge of these shape factors enables the determination 
of an<i appearing in Eqs. (2), (3) and (8), by the 

use of standard expressions for radiant energy exchange between 
coaxial cylinders with gray inner surfaces \ 22 j . The values 
of 3? g calculated in this way are : 

P 1-2 = 0.802', and P = 0.895. 

The values of trsnsnittivities , absorptivities and 
amiss ivities of glass and metal surface (with a selectively 

absorbing black coating), are listed below in Table 2-1, 

TABLE 2-1 

Transmittivities, Absorptivities & Emissivitiea ;{7] 



Glass 


Single 

0.88 

0.05 0.06 

0.94 

Double 

0.81 


Aw 


Metal Tube 


0.0 1.0 0.91 0,09 


2.3 METHOD OP SOLUTION 

The coupled and non linear differential equations in- 
volved in the determination of T^, Tg, T g and T q , preclude the 
possibility of obtaining analytical solutions. The simplest mode 
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of attack is therefore to. obtain numerical solutions for para- 
metric -values of the tube diameter, D , the mass flow rate, m 
and ihe tube length. Using "the data as specified above along 
with the specified diameters which are 11 cms, 15 cm and 17»8 cm 
for metallic tube, inner glass shield and outer glass shield 
respectively, and the mass flow rate of the liquid through the 
tube, an iterative, step-by-step calculation may be carried out. 
The steps involved in the computation are the following : 

a) Choosing a sufficiently small x, say 1 cm., and assuming the 
metal surface at the inlet, equal to that of the entering 
water, the simultaneous equftions for T^ and T 0 at the inlet 
are determined numerically. The technique used in the cal- 
culations is Newton's method of tangents, assuming that the 
initial guesses are sufficiently close to the roots [ 23^ . 

b) Using Eqs.(2)-(8) as well as T 1 and T g as determined above, 

a first estimate of T 4 is obtained. In this calcula- 

s 

tion it has been assumed that x is approximately equal to 

x*/2, since x* is small (l3.9 x 10 ^ / a) and there is not 

— * 

likely to be much of difference between x and x* / 2. 

c) After setting I ( x* / 2 } = I (o) + T 4 x* )■ / 2 , 

s s s 

steps (a) and (b), are repeated until the newly computed 

values of T 4 x* )■ agrees with that computed early closely. 
s 

The values of T^ , and T 2 corresponding to T g ( x* / 2 ) and T g 
4 x* as well as the liquid temperature, T^ 4 x* ^ are de- 
termined by using the appropriate equations. 
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d) The efficiency of the solar water heater in the interval 
0 to x is 

" V ( = ^ c p ( T^ 4 x* } - Te ) / (D 2 q x) * (l6) 

e) Now, the origin may be translated to the point x = x*/2, so 
that the new starting values for the computations are 

T s ^ = T s ^ x * /2 * * (°) = ^ ( x*/ 2 ) a »d T 2 (0) = 

T^ ( x*/2 ). The calculations at steps (b), (c) and (d) are 
repeated, until the temperatures along the tube surface as 
well as the mixed mean fluid temperatures, and the efficien- 
cies at all values of x, are obtained. 

The computer programme to evaluate the temperatures 
using the above scheme is listed in Appendix I. 

The computations may be stopped, once the end of the 
tube is reached which is 1.20 m in the presaat set up. However, 
if the purpose of the calculation is to determine the solar heater 
length needed for maximum energy collection when the tube diameter 
and the mass flow rate are specified, it is necessary to continue 
the calculations until the newly computed fluid temperatures 

{ x* } does not differ appreciably from that one step computed 
earlier. For the tube diameter of 18 cms and mass flow rates of 
0.5, 1.5, 2.5, 3.5, 4.5, 5.5 and 6.5 kg/min. ' the optimum length 
has been found £24 J . For mass flow rate of 0.5 kg/min. ‘ the 
optimum length comes out to be 6 m, the difference between T at 
6 m and that at 7 m being less than 2$. 


* Please see page 24 



In the present work the computations have been done 
for mass flow rates of 0.5, 1.5,a2.5 and 3.5 kg/min. The 
comparison of results between those computed and found experi- 
mentally are shown in Pig. 4.6. 


* Hie expreaeie* i> mt valid f#r utwa! ee w rce tiea mm* the efficient] 
•f teeter la the uiwtl eMmtim >odt lit 

\ - . « f Cf f - f t ) / 1 

wfcmwm at m et# imt&tom-w Int isjfodi 'hkm .M 4k.m m 

Waa m » W Mi Mr wwmmrn mmmmmmw w» nn»MW m mm Ma^r m 

I* m fcaperaiure *f water et ttw wad mt day 

fj * Teaperature mt wter is the aoniag 

9 - ?«tol rediatica (direct, diffuse and pend reflects) 
incident « tte heater dniag the day. 




CHAPTER III 


EXPERIMENT AI SET UP 

In order to check: the theory and also to study cylindrical 
solar water heaters more thoroughly, a heater has been fabricated 
and experimented upon. The line diagram of the set up is shown an 
3.1. The main parts of the set up will be described under 
two headings, i) Apparatus and ii) Instrumentation. 

3.1 APPARATUS 

i) G.I. Tube 

A tube 2.5 m long, 10.15 cm. inner diameter, and outer- 
diameter 11 cm., made of Galvanised Iron is used for heating the 
water. In the theoretical analysis of the heater it has been 
assumed that the fLow of water in the tube is fully developed. 

Hence a sufficiently long entry length must be provided at the 
tube inlet to ensure a fully developed flow. 

The hydro dynamical entry length can be calculated by 
using the equation s {j25j 

x = 0.03 Re D 

vdtiere x = entry length, Re = Reynolds number and D = diameter 
of the tube. The mass flow rates for which the heater has been 
experimented upon vary from 0.5 kg/min to 3.5 kg/min . . Eor 3.5 
kg/min. mass flow rate, the entry length oomes to 3.2 m. Since 
such a long tube was not available the longest tube available, 
inhich is 2.5 m lcng has been used. 
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The heated section of the tube is 1.2 m. Thus the 
remaining 1.3 m length of the tube acts as ihe entry length. Por 
mass flow rates of 0.5 kg/min. and 1.5 kg/min. the entry length 
as found by equation above comes to 0.45 m and 1,35 m respectively. 
Thus fully developed flow as assumed in the theory can exist only 
for flowrates of 0.5 kg/min. and 1.5 kg/min..,. 

Since this length is meant only for developing the flow, 
it has been insulated by 5 cm thick glass wool so that sun does 
not heat up this portion. By making simple calculations it has 
also bean seen "that the amount of energy input in this section, 
because of direct, diffuse and ground reflected radiation, comes 
out to be less than 1$ of that incident on the heated section. 

The portion of the tube to be heated has been painted 
black with a paint vhidi is a mixture of lamp black and black 
enamel paint. 

Water is allowed to flow from -the lower end of the tube 
which is inclined at an argle of 6° to the horizontal, as shown 
in Pig. 3.1. The flow rate has been measured with "the tsa lp . f 
of a flow meter. 

ii) Wooden Panels 

These panels are mounted on tho tube and clamped to it 
by means of 6 bolts, 3 on each one of them, as shown in Pig. 2*1. 
At the back of these paiels are the thermocouple wire connectors. 
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In order that the connectors are- not exposed to sun they have 
been shaded by white packaging material. 

iii) Glass Shields 

The two glass fields are Of.(35«f)dan«. cmjdi&Eieter 
respectively. They are made of pyrex glass with average trans- 
mittivity of 0.89. The glass shields are 2 mm thick. The air 
gap between the G.I. tube and two glass shields is 2 an and 1.5 cm 
respectively. These shields are mounted on steps made on the 
wooden panels, as shown in Pig. 2.1. Rubber gaskets have been 
placed at the rim of the panel so that contact between glass 
shield and panel is good, thus preventing any enterance of dust 
or dew. 

3.2 IN 3TRUMENTATI0 IT 

i) Taaperature Measurements 

The temperatures of water as well as the G.I. tube and 
glass shields have been measured by calibrated 24 gauge Insulated 
copper constanton thermocouples. 16 thermocouples have been 



Fig 3 2 
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are fixed on outer and inner glass shields by means of an adhesive, 
epoxy resin. 

The water inlet and outlet temperatures have been mea- 
sured by dipping the thermocouples in inlet and the outlet pipes 
as shown in Fig. 3-1. 

The ambient temperature has been measured by putting a 
thermometer, with a least count of 0. 1 °0, in shade. 

ii) Recorder 

A Leeds and North rap millivolt recorder with a least 
count of 0.05 mV has been used Id record "the temperatures. 

iii) Solar Radiation Measurements 

d 

Solar radiation has been measured by Epp]^" pyranometer 
installed by National Meteorological Laboratory Poona at the 
Civil Engineering Department of 1*1. 1. Kanpur. The charts from 
this recorder for both direct and diffuse radiation have been 
used to calculate tie amount of radiation falling on the solar 
water heater. The radiation as read from these charts is that 
incident on a horizontal surface. The radiation incident on an 
inclined surface corresponding to the inclination has been cal- 
culated therefrom by using appropriate factors as shown in 
Chapter II. 

3.3 MODES OF TEST 

The performance of the heater has been tested under two 


conditions ; 
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i) Hatural circulation 

ii) Forced circulation 
Natural Circulation 

In this mode, the pipe is filled with water in the morn- 
tog and the tap at the lower end of the tube is closed. The 
temperatures of the tube and the glass shields are recorded. At 
the end of the day, the inlet and the outlet of the tube are 
connected by a small pump (0.025 hp) through polythene pipes so 
that the water is thoroughly mixed. The final temperature of the 
water is then recorded. 

Forced Circulation 

The water is forced in the tube from below. There is 
no use of pump as the tap pressure is enough to force the water 
to the maximum flow rate of 3.5 litres/ min. Temperatures both 
for inlet and cutlet water as well as those of the tube and the 
glass shields are measured. At the end of the day the water in 
the tube is mixed by the pump in order to measure the total ther- 
mal energy gained by the system. The heater has been tested for 
four mass flow rates, 0.5, 1.5, 2.5 and 3.5 kg. /min.. 
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RESUITS AID DISOJSSIOFS 

4*1 rescuts 

The experimental phase of this study consisted of tests 
conducted during the months of May and June, 1974. luring the 
test period, a wide lange of climatic conditions were encountered; 
particularly the onset of early monsoon made it difficult to carry 
out tests on all days. Thus only the data obtained on the days 
there was uninterrupted sunshine is presented in this work. 

The heater 'was tested in two modes viz. natural convec- 
tion and forced oonveotion. Table 4.1 shows the experimental 
data and results for natural convection. 

Table 4.1 : Experimental Data and Results for Test in 
Natural Convection (May 18, 1974). 


Water tempera-tare] 

°c 1 

[ Time i 

^ Insolation 
angle) 

(koal) i Hea-t 

v J (gained (ciency 

Initial | 


IFpom I To \ 
1A.M. 1 P.M.l 

jDirectjDiffusei&round jtby wa-terj % 
f X ^reflected? (kcal\ $ 

31.5 

55.3 

7.20 6.10 

8° 1040 487.15 

172.5 535 31.5 


Figure 4.1 shows the plot of the experimental data for 
natural convection. The direct and diffuse solar radiation fall- 
ing on the tilted collector are plotted, together with -the plate 
temperature and the ambient temperature. The maximum plate tempe- 
rature of 63.5 °C was obtained at 3.00 p.m. after which it started 




Insolation,Kcat /m 



IUU 


Time of day, hours. 

• o Direct normal radiation on c o Hoc tor surface 
4 Diffuse sky 
• Plate temperature . x 
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to decrease since the heat losses from the plate to the surround- 
ings were more than the input by solar radiation. 

The efficiency obtained in this case is rather small 
being 31*5$. A possible reason is that since the angle of tilt 
of the heater to the horizontal is small (8°), convection 
currents do not develop readily. Also since the upper surface 
of tile tube is heated there is no radial convection. More vigo- 
rous convection could have been obtained if the angle of tilt 
had been around 30° , but then this would have been at the ex- 
pense of less solar input to the heater, because then the solar 
input would not have been normal to the surface of the heater. 

Since the efficiency in natural convection is small 
a better indication of the effectiveness of the heater would be 
to test it under forced convection. This does not mean that 
the efficiency of such a type of heater will be small in natural 
convection, but it simply means that under the conditions of 
present experiment when the inclination angle is 8°, it will be 
small. Table 4.2 tabulates the results for test in continuous 
flow for different mass flow rates. 

Figures 4.2 - 4.5 show the plot of experimental results 
for forced convection. Together with the total radiation falling 
on the heater, the outer glass temperature and the difference 
between the water temperatures at the outlet and the inlet are 


plotted 
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TABLE 4.2 : Experimental Bata and Results far test in Continuous Flow 




I I « 1 I'M . W i ' — I V-* J if I 'TV '.V* I ■ 

o Direct insolation on the surface of collector cp =6° 

I%p0 • Diffuse insolation » ” ” = 6°: 

x Outer glass temperature 
D Ambient temperature 

a Difference between water temperature at outlet & inlet 

, ' : ' ' -cj -v A " - , 

FIG. 4-2 Experimental results for test at Kanpur, 29 May, 1974 for 

„ . mass flow rate of 0 ^ 
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In all these figures it is seen that the difference 
between the outlet and the inlet temperatures of water at the 
initial stage of the experiment is negative and that the time in 
which this difference remains negative decreases as the mass 
flow rate increases. A possible expl&ination for this can be 
given as follows. In the morning the heater is filled with 
water at a particular temperature. The inlet water is fed in 
from a tap which is connected by a long system of G.I. pipes 
to an overhead tank. Thus the inlet water is warmer than that 
in the heater pipe. In the early stages of the experiment, 
this hot incoming water gets mixed with the colder water inside 
the heater and this becomes cooler. Also in the morning Hie 
amount of sun’s radiation is not enough to provide as much energy 
as it has lost in the insulated portion. For the lowest mass 
flow rate of 0.5 kg/min. the time during which the outlet and 
inlet water temperature difference remains negative is about 
3 hour. As the mass flow rate increases the difference between 
temperatures of the incoming water and that in insulated section 
goes on decreasing. Also the time of contact of the incoming 
water w'ith that in the insulated section reduces. Therefore the 
time during which the difference between the outlet and inlet 
water temperature remains negative, decreases. For mass flow 
rate of 3.5 kg/min. this time reduces to less than an hour. 

From Table 4.2 it is also seen that the heater effi- 
ciency increases with the increase in mass flow rate. As the 
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mass flow rate is increased the heat transfer coefficient also increases 
and thus there is better transfer of heat from the tube to the fluid* Also 
at higher mass flow rates the temperatures of the plate are lower than those 
at lower mass flow rates and thus there are smaller losses from the plate 
to outside. A combination of these two effects loads to the observed rise 
in efficiency. If a second order regression curve is drawn through the data 
points of efficiency vs. mass flow rate then all the points lie within + 5 % 
of this curve. At lower mass flow rates of 0.5 and 1.5 kg/min the trend 
shows an increase in efficiency with increase in mass flow rates because of 
reasons stated above. For higher mass flow rates the trend is difficult to 
predict since besides the scattering of points due to experimental error 
there is the problem of flow being not fully developed because of insuffi- 
cient entry length. 

Figure 4.4 shows the experimental data for a mass flow rate of 
150 kg/h. There is a sudden dip in the incident radiation of the sun at 
about 1 p.m. This was because of clouds. This effects the water temperature 
difference whi ch also dips at 1,30 p.m. 

Finally the experimental results are compared with those obtained 
theoretically. Figure 4,6 shows the comparison for efficiency, plate 
temperature and difference between water temperature at the outlet and 
inlet, for varying mass flow rates of 0.5 kg/min. to 3.5 kg/min. 

f : ■ i' ! - ‘ -I • ■ '• . . .. , 

Table 4.3 lists the -various quantities, both: experimental 
and theoretical. The theory compares reasonably well with the 
experimental results obtained. The difference between the theore- 
tical and experimental results for mass flow rates of 30 kg/h to 
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150 kg/h lies within 1 "]%. The theoretical efficiency is higher 
since we have assumed in the analysis that the absorber plate 
is coated with a selectively absorbing paint which has absorp- 
tivity of unity and very low enissivity. In the actual test 
the emissivity of the laup black coating is also nearly equal 
to unity. Hence one may expect higher losses due to radiation 
and a reduced efficiency. Also we have assumed in theoretical 
calculations that the tube temperature is uniform at any given 
section whereas this cannot possibly occur in practise owing 
to non-uniform irradiation over the tube periphery, and the 
relatively small thermal conductivity of the tube. 


Table 4.3 J Comparison of Experimental and Theoretical 
Results of Heater under Ibrced Convection 
Mode. 


Mass flow \ Efficiency $ 

rate kg/h 1 

^Theoretical j( Exptl 


30 

74.5 

67.0 

54.0 

49.5 

4.12 

3.4 

90 

78.4 

70.3 

46.5 

44.7 

1.58 

1.29 

150 

78.35 

76.4 

45.3 

44.0 

1.05 

0.94 

210 

78.3 

62.0 

44.5 

43.5 

0.64 

0.51 


j[piate Temp. °C ^Diff. of' outlet and " 

| jjinlet water temp. °C/h 

fTheoreticallExptlf Theoretical K Exptl. 


It is worthwhile to comment here on the efficiencies 
reported by Vincze j 12 j . He has reported efficiencies as high 
as 216$. There are two probable reasons for such absurdly high 
efficiencies. One is that Vincze has assumed the radiation falling 
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The theory compares reasonably well with the experi- 
mental results. The smallest difference of 2.8 % is found for 
mass flow rate of 210 kg/h whereas for mass flow rates of 30 
and 90 kg/h the differences lie in the range of 10 - 11$. 

Still from purely heat transfer point of view there is ample 
scope of improving the theory, and it may be taken up as the 
next stage of investigation. The variation of solar radiation 
with time, together with the circumferential temperature varia- 
tion may be incorporated in the refined theory. 

The present study is also a first step towards making 
an efficient solar collector device which can be used as a gene- 
rator of a vapor absorption refrigeration system. The tempera- 
tures obtained in the present study are not high and in order 
to obtain higher temperatures a parabolic concentrator can be 
put below the cylindrical heater. The cylindrical heater has 
one more advantage over the flat - plate collectors, if used 
for purposes of generator, and that is -that the cylinder can 
withstand higher pressures easily while in flat plate collectors, 
bulging of the absorber plate and consequent leakage may set in, 
if high pressures are used. 
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APPENDIX I 


COMPUTER PROGRAM 




w 31 AUGUST r> 

IB FTC MAIN 

dimension VM( 5»,am(9 ) 

COMMON X,TS»TGT,TGj» r 
PEAL NETA1120C) 

READ, I V M ( I ) , 4 M { I ) ,1=1 
R S = 0. 5 

31 X = C13.9E-5)/RS 
READ, TS, TGI, TGJ 


SHG, TGIHC, AF,AQ,TS 


AQ = 4 . 

DO 20 1 = 1 , 5 

AP=AP+{ gXP ( ( “ VM ( I ) )*X )■" . )/ { VM( I )*AM ( I ) ) 
AQ=AQ+l./( VM( I ) #AM ( I } ) 

TE = TS 
TM=TS 

CALL ALGEQfTS ,TGI,TGJ J 

TSN = TS 

TGIN=TGI 

TGJN=TGJ 

TSI*TS + G(X,TS,TGIM{AP+AQ) 

CALL ALGEQCTS I, TGI, TGJ) 

TSH=0.5*(TSI»TSN J+TSM 

TGIH=0.5*<TGI~TGIN)+TGIN 

TG JH=0.5#( TGJ--TGJN )+TG JN 

TGIH0=TGIH 

TSHO=TSH 

TG JHO = TG JH 

CALL DIFEQ 

TSH=0.5*<TSAN*»TSNHTSN 
CALL ALGEQ(TSH,TGI »TG J) 

T6lH0=TGI 

T6JH0=TGJ 

T$HO=TSH 

*-v ■ • ‘ - . 

KRV* ; 

CALL ALGfiQ(TSH,TGI,TG J) 

TGIH=TGI 
TG JH=TG J 

GO TO 1,3 r: 

TGI=TGIHO 

TGJ«TGJH0 


ALGEQ (TSAN, TGI ,TG J) 


’GJHO^TGsl 









